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Clustering, Bias and the Accretion Mode of X-ray selected AGN 
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ABSTRACT 

We present the spatial clustering properties of 1466 X-ray selected AGN compiled from the 
Chandra CDF-N, CDF-S, eCDF-S, COSMOS and AEGIS fields in the 0.5-8 keV band. The 
X-ray sources span the redshift interval < z < 3 and have a median value of f = 0.976. 
We employ the projected two-point correlation function to infer the spatial clustering and find 
a clustering length of ro = 7.2 ± 0.6/i _1 Mpc and a slope of 7 = 1.48 ± 0.12, which corre- 
sponds to a bias of 6(f) = 2.26 ± 0.16. Using two different halo bias models, we consistently 
estimate an average dark-matter host halo mass of A/h — 1.3(±0.3) x 10 13 /i~ 1 Af©. The X- 
ray AGN bias and the corresponding dark-matter host halo mass, are significantly higher than 
the corresponding values of optically selected AGN (at the same redshifts). The redshift evo- 
lution of the X-ray selected AGN bias indicates, in agreement with other recent studies, that 
a unique dark-matter halo mass does not fit well the bias at all the different redshifts probed. 
Furthermore, we investigate if there is a dependence of the clustering strength on X-ray lu- 
minosity. To this end we consider only 650 sources around z ~ 1 and we apply a procedure 
to disentangle the dependence of clustering on redshift. We find indications for a positive de- 
pendence of the clustering length on X-ray luminosity, in the sense that the more luminous 
sources have a larger clustering length and hence a higher dark-matter halo mass. In detail we 
find for an average luminosity difference of S log 10 L x ~ 1 a halo mass difference of a factor 
of ~3. 

These findings appear to be consistent with a galaxy-formation model where the gas 
accreted onto the supermassive black hole in intermediate luminosity AGN comes mostly 
from the hot-halo atmosphere around the host galaxy. 

Key words: galaxies: active : clustering- X-rays: galaxies 



1 INTRODUCTION 

One of the most remarkable astonomical findings of the last 
decade is the discovery of the scaling relations between the 
mass of supermassive black holes (BHs) and the properties of 
the large-scale environment of their host galaxies. In particu- 
lar, observations suggest that there is a tight correlation be- 
tween the BH mass and the bulge velocity dispersion, stel- 
lar mass and lumi nosity of the spheroidal component of th e 
host galaxy (e.g., [Fg rrarese & Ford 2005 Mago rrian et al. 19981 ). 
jHaring & Rix 2004Gultekin et al. 20091) . These observational re- 
lations indicate that the cosmic growth of BH mass is strongly cou- 
pled to the evolution of the host spheroid. However, the physical 
mechanism that shapes these relations is still unknown. 

In m ost semi-analytic models of galaxy formatio n 
(e.g., ICroton et al. 200dBower et al. 200dLago's et al. 20081) . 



ISomerville et al. (20081) it is assumed that the active galactic 
nucleus (AGN) associated with accreting black hole regu- 
lates the formation of new stars. The mechanisms adopted 



by semi-analytical models for triggering AGN activity in- 
clude major galaxy mergers for the most luminous AGN 
dDi Matteo et al. 200 fepkins et al. 200dMarulli et al. 20091) . 
while it is possible that in the lowest luminosity AGN 
regime secular disk instabilities or minor interactions pla y 
the key role (Hopkins & Hernquist 200(feournaud et al. 20 TTh . 
These different AGN fueling modes make diverse predic- 
tions for the envi ronment of the galaxies that host AGN 
JShankar et al. 20091) . For example, in the major-merger sce- 
nario, only weak luminosity dependence on clustering is expected 
faopkins et al. 200 jLidz et al. 200dBonoli et al. 2009n . Hence, 
observational studies of AGN clustering and its dependence on 
luminosity can place valuable constraints on the AGN fueling 
modes and consequently on the AGN-galaxy co-evolution models. 

The clustering of AGN has been studied with excel- 
lent statistics mainly in the optical and particularly in large 
area surveys, such as the 2QZ (2dF Q SO Redshift Survey 
ICroom et al. 200 jPorciani & Norberg 20061) and the SDSS (Sloan 
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Digital Sky Survev. lLi et al. 200dShen et al. 200 jRoss et al. 20091) . 
These surveys found no evidence for a strong 
dependence of clustering on luminosity (e.g., 
ICroom et al. 200 jMvers et al. 200rjda Angela et al. 20081) al- 



though [Sheri_eXaL_(2Q09|) detect an excess of clustering for their 
10% brightest quasars. However, optical QSO may represent only 
the tip of the iceberg of the AGN population. Ver y deep X-ray sur - 
veys find a surface density of about 10,000 deg -2 l lXue et al. 20 111 ), 
which is about two orde rs of magnitude higher than that found 
in optical QSO surveys dWolf et al. 20030 . Several studies have 
explored the angular clustering of AGN in X-ray wavelen gths 
Jvikhlinin & Forman 1995kkvlas et al. 200(lYang et al. 2003h. 
jBasilakos et al. 2004Basilakos et al. 20051) . iGandhi et al. (20061) 



Puccetti et al. (20061) 
Plioni s et al. (20081) . 



lEbrero et al. (2009) . lElviv et al. (20121) . 2 AGN CATALOGS 



Carrera et al. (20071) . iMivaii et al. (2007) 



our findings for the AGN correlation function as estimated by the 
joint and individual X-ray samples. In Section 5, we calculate the 
bias of AGN and its redshift evolution, and estimate the average 
mass of the DM halos hosting X-ray AGN. In Section 6, we in- 
vestigate the dependence of the halo mass on X-ray luminosity and 
compare our findings with theoretical galaxy formation models for 
the evolution of accreting BHs. Finally, we present in Section 7 
our conclusions. For comparison reasons with previous works we 
adopt, unless otherwise stated, a flat cosmology with a present- 
day matter density parameter Q. m — 0.3, a cosmological constant 
£Ia = 0.7 and baryon density Qt = 0.04. The Hubble constant is 
expressed in units of h as Ho = 100 h km s _1 Mpc . 



These studies measure t he projected a ngular clustering and then 
via Limber's equation l lPeebles 1980h derive the corresponding 
spatial clustering length. However, in this method an a priori 
knowledge of the redshift distribution d N/dz is needed and thus, 
the uncertainties may be appreciable. IPlionis et al. (20081) first 
reported strong indications for a luminosity dependent clustering 
of X-ray AGN in the CDF fields for both soft and hard bands in the 
sense that X-ray luminous AGN lie in more massive dark-matter 
(DM) halos compared to the less luminous ones. 

Recently, several studies have attempted to measure 
the spatial correlation function of X-ray selected AGN by 
employing sp ectroscopic redshifts to estimate their correc t 
distances <Mullis et al. 2004Gilli et al. 200 feme et al. 20060 . 



dGilli et al. 2009hickox et al. 200 ^Coil et al. 20090. 
IKrumpe et al. (20ld) JCappelluti et al. 20ldMivaii et al. 20110 , 
dStarikova et al. 201 llAllevato et al. 201 lh . These studies show 
that X-ray AGN are typically hosted in DM halos with mass of 
the order of 12.5 < log lo (M h /[/i _1 M ]) < 13.5, at a redshift 
z ~ 1. On the issue of luminosity de pendent clustering the 
results are contentious. lYang et al. (20061) did not detect strong 
correlation between X-r ay luminosity and the clustering amplitude 
for Chandra sources, fcilli et al. (20091) performed a similar 
investigation in the XMM-COSMOS field dividing the sample 
below and above L x (0.5 — lOkeV) = 10 44 erg s -1 and did no t 
find any luminosity dependence. Similarly. IStarikova et al. (201 lh . 
using 1282 Chandra/Bootes sources, did not find a significant 
dependence of clustering on luminosity. 

On the other hand, ICoil et al. (20091) using AGN in the red- 
shift range z — 0.7—1.4 from the AEGIS field found a weak 
evidence for a luminosity depende nt clustering, but not a t a statisti- 
cally significant level. In addition. IKrumpe et al. (20100 using low 
redshift (z ~ 0.25) AGN, selected from the ROSAT all-sky survey 
and cross matched with the SDSS, found that L x > 10 44 erg s _1 
sources are clustered m ore strongly than lower luminosity sources. 
In another recent study, ICappelluti et al. (20 id) using a sample of 
199 Swift-BAT sources in the 15—55 keV band, found a marginally 
significant luminosity dependent clustering. 

In this paper we use a sample of 1466 sources selected in the 
0.5 — 8 keV band from a variety of deep Chandra X-ray surveys, to 
estimate the spatial correlation function and typical DM host halo 
mass of X-ray selected AGN. Our aim is to investigate whether 
there is a luminosity dependence on clustering, but also to obtain 
insight into the black-hole fueling mechanisms, using the largest 
X-ray selected AGN sample used so far for this scope. The paper 
is organized as follows. In Section 2, we briefly discuss the X-ray 
surveys we consider in this paper and in Section 3, we describe the 
methodology of our clustering analysis. In Section 4, we present 



We make use of X-ray data coming from the five deepest X-ray 
surveys, namely the Chandra Deep Field South and North (CDF- 
S and CDF-N), the AEGIS, the extended Chandra Deep Field 
South (ECDF-S) and the COSMOS survey. These cover a vari- 
ety of exposures (from 4 Ms down to 40 ksec) and surveyed area 
and thus, cover extensively the luminosity-redshift space. More- 
over, these fields contain excellent quality spectroscopic observa- 
tions and therefore good quality redshifts, which are essential in 
this project. Below we present briefly the main characteristics of 
the X-ray surveys used in this work. 

2.1 CHANDRA DEEP FIELD NORTH 

The deep pencil beam CDF-N survey covers an area of 
448 arcmin 2 , is centered at a = 12 h 36 m 49 s ,<5 = +62°12'58" 
and consists of 20 individual ACIS-I (Advanced CCD Imag- 
ing Spectrometer) pointings. The combined observations pro- 
vide the deepest X-ray sample currently available together with 
the CDF-S. Here, we use the X-ray 2 Ms source catalogue of 
Alexander et al. (20031) . with a sensitivity of ~ 10 erg cm 
s _1 , which consists of 503 sources in the 0.5—8 keV band. Spectro- 
scopic redshifts were used for 243 X-ray s ources in the redshift in- 
terval z = — 3 from lTrouille et aL (20081 and references therein). 

2.2 CHANDRA DEEP FIELD SOUTH 

The deep pencil beam CDF-S survey covers an 
area of 436 arcmin 2 and the average aim point is 
a = 03 h 32 m 28 s 8, S = -27°48'23"(J2 000). The an alysis of 
all 23 observations is presented in l lLuo et al. 20081). We use 
the 2Ms X-ray source catalogue of iLuo et al. (20081) . which 
consists of 462 X-ray sources. Spectroscopic redshifts were used 
for 219 X-ray s ources in the redshift interval z = — 3 from 
iLuoet al. (20100 . 

2.3 AEGIS 

The ultra deep field survey comprises of pointings at 8 sepa- 
rate positions, each with nominal exposure of 200 ks, cover- 
ing a total area of approximately 0.67 deg 2 centered at a = 
U h 17 m ,S = +52°30' in a strip of length 2 deg. The flux limit 
of the survey is S ~ 10 -16 erg s -1 cm~ 2 in the f ull band. We 
use the X-ray source catalogue of lLaird et al. (20091) , which con- 
sists of 1325 sources. Spectroscopic redshifts were used for 392 
X-ray sources in the redshift interval z = — 3 from DEEP2 
dPavis et al. 200llDavis et al. 200jCoil et al. 20091) . 
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2.4 COSMOS 

The Chandra COSMOS Survey covers the central 0.5 deg 2 area of 
the COSMOS field with an effective exposure of ~ 160ks and the 
rest of the field with an effective exposure of ~ 80ks. The limiting 
source detection depths are 5.7 x 10~ 16 erg s -1 cm " 2 in the full 
band. We use the source catalog of lElvis et al. (20091) . which con- 
sists of 1761 sourc es. Spectroscopic redshifts were used for 417 
X-ray sources from lBrusa et al. (20101) . 



2.5 EXTENDED CHANDRA DEEP FIELD SOUTH 

The ECDF-S survey consists of 4 Chandra 250 ks ACIS-I point- 
ings covering ~ 0.3 deg 2 and surroundin g the original CDF - 
S. S ource detection ha s been performed by iLehmer et al. (20051) 
and IVirani et al. _(2006l) . Here we use the source detection of 
ILehmer et al. (20051) in which 762 sources have been detected. The 
flux limit of ECDF-S is S ~ 10~ 16 erg s _1 cm -2 . Spectroscopic 
redsh ifts were used for 288 s ources in the redshift interval z = 0—3 
from Isilverman et al. (20ld) . When combining all fields, we ex- 
clude sources that are detected both in ECDF-S and CDF-S and 
keep only those detected in the ECDF-S. 



3 METHODOLOGY 

3.1 Theoretical Considerations 

The main statistic used to measure the clustering of extragalactic 
sources is the two-point correlation function £(r). £(r) describes 
the excess probability over random of finding a pair with one ob- 
ject in an elemental volume dVi and the seco nd in the elem en- 
tal volume (IV2, separated by a distance r (e.g.. |Peebles 1980h . Its 
mathematical description is given by dP = (n) 2 [l + £(s)]dVidV2, 
where (n) is the mean space density of the sources under study 
and s = cz. When measuring £ directly from redshift catalogues 
of sources, we include the distorting effect of peculiar velocities, 
since the true distance of a source is r — (s — v p • r)/ Ho, where 
v p • r is the component of the peculiar velocity of the source 
along the line of sight. In order to avoid such effects one can ei- 
ther measure the angular clustering, which is not hampered by the 
effects of ^-distortions, and then, under some assumptions, infer 
the s patial correlat ion function through the Limber's integral equa- 
tion ^Limber 19531) . Alternatively one can use the redshift informa- 
tion t o measure the so-called projected correlation function, w p (r p ) 
(e.g jDavis & Peebles 19831) and then infer the spatial clustering. 

To this end, one deconvolves the redshift-based distance of 
a source, s, in two components, one parallel (n) and one perpen- 
dicular (r p ) to the line of sight, i.e., s = (r 2 + 7r 2 ) 1//2 , and thus 
the redshift-space correlation function can be written as f(s) = 
£(r p , 7r). Since redshift space distortions affect only the it com- 
ponent, one can estimate the free of 2-space distortions projected 
correlation function, w p (r p ), by integrating £(r p , n) along n: 



w p (r p ) = 2 / C( r pj 7r)d7r. 



(1) 



Once we estimate the projected correlation function, w p (r p ), 
we can recover the re al space correlation fun ction, since the two are 
related according to jPavis & Peebles 19831) 



w p (r p ) = 2 J™ C(\/r 2 +^ 2 )d^ = 2 £ j| 



r£(r)dr 



Modelling £ (r) as a power law, £ (r) = (r/ro)~ 7 one obtains, 



(2) 



with 



(3) 



(4) 



where T is the usual gamma function. 

However, it should be noted that although Eq. ([3} strictly holds 
for 7r max = 00, practically we always impose a cutoff 7r max (for 
reasons discussed in the next subsection). This introduces an un- 
derestimation of the underlying correlation function, which is an 
increasing function of separation r p . For a power law correlation 
function this u nderestimation is eas ily inferred from Eq. $2^ and is 
given by (e.g.. lStarikova et al. 201 ll) 



Cj(r p ) 



f 7T n 

Jo 



J °°(r 2 -7r 2 )-7/2d7T 



(5) 



Thus, by taking into account the above statistical correction, and 
under the assumption of the power-law correlation function, one 
can recover the corrected spatial correlation function, £(r p ), from 
the fit to the measured w p (r p ) according to (which provides also 
the value of 7): 



£(r,) = 



u p (r p ) 



-A7 C*y ( Vp ) Vp 



(6) 



However, at large separations the correction factor increasingly 
dominates over the signal and thus it constitutes the correction pro- 
cedure unreliable. Alternatively, as can be easily shown using Eqs. 
® and (6), a crude estimate of the corrected spatial correlation 
length can be provided by: 

r o ,c = r C 1 {r r lh , (7) 
where ro and 7 are derived from fitting the data to Eq. {3). 



3.2 Correlation Function Estimator 

As a first step we calculate (,(r p ,ir) using the estimator (cf., 
iKerscher et al. 2000l) 



£(r p ,7r) = 



N R DD(r p ,n) 
N D DR(r p ,n) ' 



(8) 



where Nd and Nr are the number of data and random sources, 
respectively, while DD(r p ,n) and DR(r p . n) are the number of 
data-data and data-random pairs, respectively. We then estimate 
the redshift-space correlation function, £(s), in the range s = 
0.16 — 40 hT 1 Mpc and the projected correlation function, w p (r p ), 
in the separation range r p = 0.16 — 20 h^ 1 Mpc. Note that large 
separations in the n direction add mostly noise to the above estima- 
tor and therefore the integration is truncated for separations larger 
than 7T ma x. The choice of 7r ma x is a compromise between having 
an optimal signal to noise ratio for £ and reducing the excess noise 
from high n separations. The majority of studies in the literature 
usually assume 7r max 6 [5, 30] h -1 Mpc. 

The correlation function uncertainty is estimated according to 



V3(l + w v )/\ 



(9) 



which correspo nds to that expected by the bootstrap technique 
dMo et al. 1992b . In this work, we bin the source pairs in logarith- 
mic intervals of <51og 10 (r p , n) = 0.15 and <51og 10 (s) = 0.17 for 
the w p (r p ) and f (s) correlation functions respectively. Finally, we 
use a x 2 minimization procedure between data and the power-law 
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OJ ] 10 0.1 1 10 

r p (h- 1 Mpc) s (h-1 Mpc) 

Figure 1. Comparison of the correlation function results based on the G05 
random construction method (filled symbols) with that of the standard 'sen- 
sitivity map' method (empty symbols), applied on the AEGIS full band 
(0.5-8 keV) catalog. Uncertainties correspond to 1<t Poisson errors. The 
best fit power laws are shown as black lines for the G05 method and as 
dashed red lines for the sensitivity map method. In the inset panels we show 
the likelihood lcr, 2a, 3<r contours in the 2-parameter solution space. Right 
Panel: Redshift space correlation function, Left Panel: Projected cor- 
relation function, w v (r p ). 



model for either type of the correlation function to derive the best 
fit ro and 7 parameters. We carefully choose the range of separa- 
tions in order to obtain the best power-law fit to the data and we 
impose a lower separation limit of r v ~ 1.5 hT 1 Mpc to minimize 
non-linear effects. 



3.3 Random Catalogue Construction 

To estimate the spatial correlation function of a sample of sources 
one needs to construct a large mock comparison sample with a ran- 
dom spatial distribution within the survey area, which also repro- 
duces all the systematic biases that are present in the source sample 
(i.e., instrumental biases due to the Point Spread Function variation, 
vignetting, etc). Also special care has to be taken to reproduce any 
biases that enter through the optical counterpart spectroscopic ob- 
servations strategy (cf., due to the positioning of the masks within 
the field of view and of the slits within the masks, etc). 

To this end we will follow the random catalogue construction 
procedure of iGilli et al. (2005. hereafter G05), which is based on 
reshuffling only the source redshifts, smoothing the correspond- 
ing redshift distribution, while keeping the angular coordinates un- 
changed. We will test the efficiency of this method by comparing 
the outcome correlation function results with those of the standard 
method that takes into account the details of the sensitivity maps of 
each XMM pointing. We further test the two methods using N-body 
simulations for the underlying DM density field. 



3.3.2 Standard 'sensitivity map' method 

According to the standard method of producing random cata- 
logs, each simulated source is placed at a random position on 
the sky, with a flux randomly extracted from the observed source 
log TV — log S (source number-flux distribution). If the flux is above 
the value allowed by the sensitivity map at that position, the sim- 
ulated source is kept in the random sample and a random redshift 
is also assigned to it from the observed source redshift distribu- 
tion N(z) (optimally taking into account the variation of N(z) as 
a function of flux). The disadvantage of this method is that it does 
not take into account any unknown inhomogeneities and systemat- 
ics of the follow-up spectroscopic observations. 

3.3.3 Testing the efficiency of the two methods 

In order to test the efficiency of the two random catalogue construc- 
tion methods we use the AEGIS survey, consisting of 8 fields. To 
construct the s tandard method r andom catalogue we use the sensi- 
tivity maps of l lLaird et al. 20091) and we compute the projected and 
redshift-space 2-point correlation functions. The best fit correlation 
length found is ro = 4.28 ± 0.14 h" 1 Mpc for 7 = 1.8 which is 
in excellent agreement with the results based on the G05 method: 
r = 4.13 ± 0.14 h,- 1 Mpc. In Fig.[j]we compare the correlation 
functions for the two random construction methods. 



3.3.4 A further test using simulations 

We use a 500 h^ 1 Mpc cube simulation of a flat 
il m = 0.3 and as = 0-9 c osmology with 512 3 particles 
dRagone-Figueroa & Plionis 20071) , to further test the robustness 
of the G05 random sample construction method. In order to have 
a relatively wide redshift range we replicate randomly the box to 
get an effective volume of 1500 3 h~ 3 Mpc 3 . To obtain a relatively 
large but manageable source density, we selected DM halos with 
M h > 1.5 x 10 14 /i _1 M Q , which resulted in a total number of 
haloes within our large simulation volume of TV = 9 x 2312. 

We estimate the actual halo correlation function according to 
the direct method, which entails counting the number of DM halo 
pairs within spherical shells around a given halo. The correspond- 
ing number of random pairs is estimated by {nS)V, with (n) the 
mean number density of halos in the whole simulation volume and 
5V the volume of the spherical shell i.e. 4/37r<5?- 3 , with Sr being 
the width of the spherical shell. 

We also estimate the halo correlation function using the G05 
method, for which we transform the halo Cartesian coordinates into 
spherical coordinates, determining also their distance distribution. 
The clustering results of the two methods are in excellent agree- 
ment providing 7 = 1.72(±0.05) and 1.68(±0.12) as well as 
r = 13.9(±0.3) and 13.2(±0.5) h~ L Mpc, for the direct and 
G05 correlation function estimations, respectively. 



3.3.1 Random redshift reshuffling method 

In order to assign random redshifts to the mock sample, the source 
redshift distribution is smoothed using a Gaussian kernel with a 
smoothing length of cr 2 = 0.3. This offers a compromise between 
scales that are either too small, and thus may reproduce the z-space 
clustering, or too large and thus over-smooth the observed redshift 
distribution. We verified that our results do not change significantly 
when using a smoothing length in the range er z = 0.2 — 0.4. 



4 RESULTS 

4.1 The Joint X-ray Sample 

Combining all five fields, we obtain in total 1466 X-ray sources 
with spectroscopic redshifts. The median redshift of the source 
sample is z = 0.976. In order to estimate the joint sample cor- 
relation function we generate random catalogs for each field sepa- 
rately and then combine them together. Note that we take care to 



© 0000 RAS, MNRAS 000, 000-000 



X-ray AGN Clustering 5 




■■■■■lift 



10 



20 

ax (h- 1 Mpc) 



30 



Figure 2. The resulting amplitude ro , and slope 7, of the correlation func- 
tion analysis as a function of the value of the 7r ma x cutoff. The filled circles 
correspond to the results based on w p (r p ), while the red squares on the 
results of £(r p ). The continuous red line corresponds to the crude estima- 
tion of the corrected spatial amplitude based on Eq. {7), while the dashed 
blue line correspond to the finally adopted results (which corresponds to the 
Tmax = 10 Mpc case). 



count only once the sources that are common in both the CDF-S 
and eCDF-S fields. 

We remind the reader that we fit a power-law to the measured 
correlation function over a range of separations for which such a 
model represents well the data. In order to avoid non-linear effects 
we use a lower separation limit, r p ~ 1.5 h^ 1 Mpc, for the fit of 
s. We also impose an upper separation limit (~ 10 h^ 1 Mpc) since 
we find a change above this scale in the slope of both projected and 
z-space correlation functions. 

We now investigate the sensitivity of the correlation function 
on the cutoff separation along the line of sight, 7r max , which we 
vary in the range [5, 30] Mpc. For 7r max > 30 h~ 1 Mpc the 
noise introduced by uncorrelated pairs reduces significantly the cor- 
relation function amplitude and slope. In Fig. (O we show how the 
amplitude ro and the slope 7 of the projected correlation function, 
w p (r p ) (open circles), and the corresponding real-space correlation 
function, £(r p ) (red squares), depend on 7r max . It is clear that both 
ro and 7 are relatively constant in the investigated 7r max range. Our 
best fitting values for ro and 7 are those indicated by the dashed- 
blue line (see Table 1), which corresponds to 7r max = 10 h^ 1 Mpc. 

Our main clustering results, based on the joint X-ray AGN 
sample, are shown in Fig. [3] The upper panel presents the projected 
correlation function, while the lower panel presents the redshift- 
space correlation function (circular points) as well as the inferred 
real-space correlation function, £(r p ), (magenta pentagons). The 
filled circular points are those which have been used to fit the 
power-law model correlation function, which is shown as a black 
continuous line (the dashed red line is the power-law fit for a fixed 
slope 7 = 1.8). The corresponding best fit values for the slope 7 
and the correlation length are shown in Table 1. The uncertainty 
of ro, indicated in Table 1, does not include the effects of cosmic 
variance, which can be estimated analytically, assuming a power- 
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Figure 3. Correlation functions of the joint X-ray point-source sample: The 
projected correlation function, w p (r p ) is shown in the upper panel, and the 
2-space correlation function, £(s), is shown in the lower panel. The filled 
points indicate the range over which a power-law fit was applied (black line 
corresponds to a fit with free 7, while the red line to that for 7 = 1.8). 
The real-space correlation function, £( r p) (via Eq. [6) is indicated in the 
lower panel by magenta pentagons. The inset panels show the 1, 2 and 3<r 
likelihood contours in the 2-parameter plane of the power-law fit solutions. 



law correlation function, according to: a^„ ~ ^2(7)(rp/ r) 7 , with 
J 2 (t) = 72/ [(3 - 7) (4 - 7) (6 - 7 )2 7 ] ^Peebles 1980l) . We find 
cr c „ ~ 0.3, which is smaller although of the same order as the fit- 
ting uncertainty indicated in Table 1 . 

We can single out three important results from the clustering 
analysis of the joint X-ray point-source sample: 

(i) the inferred spatial correlation length, from the power-law 
model fit to Eq. (6) or directly from Eq. (7), is ~ 20 — 25% larger 
than that estimated from the projected correlation function, a fact 
which is attributed to the correlation function underestimation im- 
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Table 1. Clustering results for the joint sample of all 5 fields (1466 
sources). The clustering length units are h^ 1 Mpc. The results correspond 
to 7r max = 10 h~ 1 Mpc, but they are very similar over the whole indicated 
Wmax range, but more so for 5< 7r max < 16 h~ 1 Mpc, as can be seen from 
Fig. E] The cosmic variance is estimated to give a further ~ 0.3 uncertainty 
in r . 



7 r r (7 = 1.8) 

w p (r p ) 1.84±0.07 5.2±0.5 5.2 ± 0.5 

£0) 1.49±0.20 6.6±1.0 6.0 ± 0.8 

1.48±0.12 7.2±0.6 6.5 ± 0.4 



1000 e 




Q ^ I I III | I | I I I I I | 

1 10 

r p ( hr 1 Mpc) 

Figure 4. Correlation function of the individual X-ray surveys: CDF-N 
(green pentagons), CDF-S (red filled squares), eCDF-S (black filled trian- 
gles), COSMOS (magenta empty squares), AEGIS (blue empty triangles) 
for redshift interval < z < 3). Upper Panel: Projected correlation func- 
tion. Lower Panel: The resulting spatial correlation function. 

posed by the necessary cutoff along the n direction in the projected 
correlation function measure, 

(ii) it appears that redshift-space distortions do not significantly 
affect the z-space correlation function, since £(s) ~ £(>>), and 

(iii) the real-space correlation function, at a median redshift of 
0.976, has a slope 7 ~ 1.5 and a correlation length ro — 7.2 h^ 1 
Mpc. 



4.2 The Individual Field Correlation functions 

We now perform the clustering analysis from the previous section 
in each individual field and compare our results to those from the 
original studies. To this end we use the redshift interval < z < 3 



and sources with luminosities L x ^ 10 41 erg s _1 , in order to 
avoid the gross contamination of our X-ray AGN sample by normal 
galaxies. We estimate the projected and z-space correlation func- 
tions in ten logarithmic intervals of width 5 log 10 (r p , n) = 0.175 
and S log 10 (s) = 0.2, respectively. With the exception of the CDF- 
N field (see below), the power-law model fit was performed for 
r p > 1.5 h~ 1 Mpc in order to avoid the non-linear contributions to 
the correlation function. 

In Fig.|4]we show the projected (upper panel) and the inferred 
via Eq. (6) spatial (lower panel) correlation function of the individ- 
ual fields, color coded, while in Table 2 we summarize their best 
power-law fit parameters 7 and ro. It is evident that the individual 
determinations of the correlation function are hampered by cosmic 
variance. It is also interesting to note that the results of the largest 
survey (COSMOS with ~ 1 deg 2 area) are consistent with those of 
the Joint sample (see Table 1). 

We find small differences with respect to the results of the 
original studies, which are mostly due to the slightly different 
sample definitions and different choices of 7r max . For example, 
for the CDF-S field we use the 2 Msec catalog rather than the 
1 Msec catalog used by iGilli et al. (20051), while fo r the CDF-N 
dGilli et al. 20051) and for the AEGIS dCoil et al. 20091) samples, the 
corresponding authors used sources with L x ^ 10 42 erg s _1 . When 
using the latter luminosity limit and the original value of 7r max , we 
recover the original clustering resultfl The only field for which it 
was not possible to fit a power-law correlation function in the linear 
scales is the CDF-N field, which presents a knee at r p ~ 3.5 h~ l 
Mpc. Therefore, we have performed a power-law fit in the two 
ranges where such is exhibited (i.e., for r p <J 3.5 and > 3.5 h^ 1 
Mpc, respectively). 

The largest clustering amplitude is observed for the CDF-S, 
a fact which has be en attributed to t he existence of a large super- 
cluster at z ~ 0.7 dGilli et al. 20031) . and which appears to affect 
also the corresponding amplitude of the eCDF-S field, the correla- 
tion function of which is estimated here for the first time. However, 
inspecting Fig. [4] and Table 2, it is interesting to note that while 
the small-scale (r v < 3.5 ft -1 Mpc) correlation function of CDF- 
N is consistent with that of both the COSMOS and AEGIS fields, 
its large-scale correlation function (r p > 3.5 h^ 1 Mpc) seems to 
be consistent with that of the CDF-S and eCDF-S fields. This sug- 
gests that either there is some unique structure affecting also the 
clustering pattern in the CDF-N (although at a lesser extent than 
in the CDF-S) or there is some unknown systematic affecting both 
CDF's. 



4.3 Luminosity Dependence of Clustering 

We now wish to investigate the possibility of an X-ray luminosit y 
dependent A GN clustering, as sugg ested by IPlionis et al. (20081) . 
Interestingly, ICappelluti et al. (20 id) , have found a dependence of 
clustering on the X-ray luminosity in the sample of local ((z) = 
0.05) X-ray selected AGN from the SWIFT/BAT survey. Similarly, 
iKrumpe et al. (201CI) have also found such indications in the X-ray 
AGN sample from the ^OSAr-All-sky-survey. 

The five X-ray surveys that we have analysed in the previous 
sections, although they cover a similar redshift range, have X-ray 
luminosity distributions which are quite different. To disentangle 

1 Due to the fact that in the luminosity dependent analysis of the next 
section we are left with a small number of sources, we are obliged to lower 
the X-ray luminosity limit to L x = 10 41 erg s — 1 . 
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Table 2. Parameters of the power-law fits to the different survey correlation functions. The results of the projected correlation function, w p (r p ), are shown 
in columns 4 and 5 and of the inferred real-space correlation function, §(r p ), in columns 6 and 7. Finally column 8 shows the correlation length of §(r p ) for 
the average value of the slope, 7 = 1.44 (which is practically equal to the slope of the joint sample, see Table 1). In all cases, except for the CDF-N, the 
power-law fits have been applied to the correlation function data, excluding the non-linear regime (< 1.5 Mpc). Since the CDF-N w p (r p ) shows a knee 
around r p ~ 3.5 h~ 1 Mpc, we present power-law fits for data below and above this scale. Correlation lengths and 7r max are in h~ 1 Mpc. 



Survey 


# 




7 


ro 


7c 


ro,c 


7=1.44 

r 0,c 


CDF-N 1 


243 


20 


1.72 ± 0.30 


4.6 ±0.7 


1.66 ±0.16 


5.1±0.6 


6.0±0.8 


CDF-N 2 


243 


20 


1.80 ±0.23 


9.9 ±3.9 


1.46 ±0.22 


13.3±4.3 


13.3±4.3 


CDF-S 


219 


25 


1.90 ±0.08 


10.6 ± 1.4 


1.72 ±0.08 


12.8±1.1 


15.1±2.1 


AEGIS 


392 


20 


1.59 ± 0.09 


4.3 ±0.6 


1.22 ±0.20 


5.5±1.0 


5.7±0.8 


COSMOS 


417 


25 


1.65 ± 0.08 


5.5 ±0.8 


1.36 ±0.15 


6.8±1.0 


7.0±0.8 


eCDF-S 


288 


10 


1.46 ± 0.07 


7.6 ±1.0 


1.23 ±0.16 


12.4±2.2 


11.1±1.8 



1 r p < 3.5 ft" 1 Mpc 

2 r p > 3.5 hT 1 Mpc 



the luminosity and redshift dependence of clustering, we construct 
a low-La; and a high-La; X-ray luminosity sub-sample for each 
survey. We also impose them to have the same redshift distribu- 
tion. This is achieved by first splitting the sources according to the 
selected luminosity limit and then matching their binned redshift 
distribution in the common area. We do so by randomly select- 
ing sources from the most abundant L x subsample so as to repro- 
duce exactly the shape of the binned redshift distribution of the less 
abundant subsample. The bin size used is Sz = 0.05. This random 
selection process is performed eight times and the final results are 
the average over these eight realizations. 

As an example of our procedure we show in Fig. \5\ the red- 
shift distribution of the whole AEGIS field (black thick line), divid- 
ing it also in its two luminosity subsamples, i.e., the low-La; with 
L x < 10 43 erg s _1 (black dense-shaded histogram) and a high- 
La; with L x ^ 10 43 erg s _1 (red sparse-shaded histogram). In the 
inset panel we also show the X-ray luminosity distribution of the 
AEGIS field, were a clear division can be observed at L x ~ 10 43 
erg s _1 . Finally, the left panel of Fig. [6] shows the redshift distribu- 
tion of the low-La; and high-La; subsamples in the common redshift 
range (0.5 < z < 1.5), while in the right panel we present one ran- 
dom realization of the matched z-distributions of the two subsam- 
ples, which are those finally used to estimate the correlation func- 
tion of the low-La; and high-La; subsamples. It is evident that our 
procedure has managed to create subsamples that have a common 
redshift distribution, therefore surpressing any redshift dependent 
effect in the comparison of the clustering results of the low and 
high-La; AGN subsamples. 

We then use the \ 2 minimization procedure to fit a power-law 
model to the measured w p (r p ) of the low-La; and high-La; sub- 
samples, typically in the projected separation range r p £ [0.16 — 
10] hT 1 Mpc. Note that it is not possible to exclude in this analy- 
sis the non-linear scales (< 1.5 h~ Mpc) due to the small number 
of sources involved and the consequent noisy correlation function. 
Furthermore, and in order to be able to compare the different re- 
sults on an equal footing, we present the inferred spatial correlation 
length, based on a power-law fit to the £(r p ) (Eq. 6) for a common 
slope, which is the mean of the low and high-La; sub-sample £(r p ) 
slopes. 

The results of the correlation function analysis for the low 
and high-La; subsamples of the individual X-ray surveys can be 
found in Table[3] In Fig.[7jwe present the inferred spatial clustering 
length, ro, c , as a function of the X-ray luminosity. There are indi- 
cations of a luminosity dependent clustering in all X-ray surveys, 
except for the CDF-S in which the correlation lengths for both low 



Table 3. Correlation function results from the X-ray luminosity separated 
subsamples for each individual X-ray survey. Note that due to the small 
samples involved we have used different bin sizes for each sample, in order 
to minimize the intrinsic scatter. Note also that the following redshift inter- 
vals were used: CDF-N 0.5 < z < 1.2; CDF-S 0.7 < z < 1.4; AEGIS 
0.5 < z < 1.5; COSMOS 0.6 < z < 1.8; eCDF-S 0.7 < z < 1.6. 
The last two columns list the slope (averaged over the two luminosity sub- 
samples) and amplitude of the power-law fit to the spatial correlation func- 
tion, £(>>), of Eq. (6). 



login L * 


# 


2 


Tmax 7 


ro 


7c 


ro,c 










CDF-N 








< 42.2 


73 


0.84 


10 


1.74 ±0.15 


3.4 ±0.4 


1.36 


4.7 ±0.7 


> 42.2 


79 


0.96 


10 


1.42 ±0.16 


4.5 ±1.1 


1.36 


7.2 ± 1.1 










CDF-S 








42.5 


53 


0.95 


10 


1.97 ±0.20 


10.0 ±2.5 


1.54 


16.5 ±2.3 


> 42.5 


44 


1.02 


10 


1.65 ±0.15 


9.3 ±2.3 


1.54 


13.6 ± 2.6 










AEGIS 








43 


90 


0.83 


20 


2.03 ±0.29 


4.6 ±1.5 


1.83 


5.1 ± 1.3 


> 43 


63 


0.82 


20 


1.90 ±0.30 


5.7±2.1 


1.83 


6.4 ± 1.8 








COSMOS 








< 43.5 


72 


1.06 


30 


2.00 ±0.13 


6.6 ±1.3 


1.85 


7.6 ± 1.5 


> 43.5 


51 


1.07 


30 


2.00 ±0.23 


10.3 ±3.0 


1.85 


11.9 ± 3.0 










eCDF-S 








< 43.2 


51 


1.04 


30 


1.91 ±0.36 


5.9 ±2.5 


1.85 


6.3 ±2.5 


> 43.2 


30 


1.02 


30 


1.98 ±0.28 


7.7 ±3.0 


1.85 


8.6 ±3.0 



and high-La; subsamples are atypically very high (prob ably due to 
the dominance of a known supercluster. lGilli et al. 20031) . Although 
the luminosity dependent clustering indications in each of the four 
remaining X-ray field are rather weak, there is an overall sugges- 
tive and systematic trend with the high La; subsamples being more 
strongly clustered than the lower La; ones. A linear least square 
fit to the data of Fig. UJ taking into account the uncertainties in 
both axes and excluding the CDF-S, provides the following depen- 
dence ro, c /(/i _1 Mpc) ~ 2.41og 10 La; — 94.7. This is indicated 
by the thick continuous line in Fig. 7. The relatively strong trend 
implied by this relation, although the individual survey luminosity- 
dependence results are not that significant, is due to the fact that 
in each field a different high and low luminosity range is sampled, 
increasing the X-ray luminosity dynamical range. 
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Figure 5. The redshift distribution of the AEGIS sources (thick black line 
histogram) divided also in the two luminosity subsamples (low-La; : black 
dense-shaded histogram, and high-La : red sparse-shaded histogram). The 
inset panel shows the X-ray luminosity distribution. The trough between 
the two apparent peaks, at L x ~ 10 43 erg s _1 , defines the division value 
between the low-La; and high-La; subsamples. 





0.5 



1.5 



Figure 6. Left Panel: Redshift distribution at the common interval z = 
0.5 — 1.5 of the Low-La; (black line) and High-La; (red shade) subsam- 
ples. Right Panel: The matched redshift distributions of the further reduced 
subsamples, according to method described in the text imposed in order to 
cancel redshift dependent effects in their correlation function comparison. 



5 BIAS AND HOST HALO MASS 

It is well established that the extragalactic sources are bi- 
ased tracers of the underlyi ng mass fluctuation field (e.g., 
iKaiser 198lBardeen et al. 1986h . The parameter that encapsulates 
this fact is the so-called linear bias factor, b, defined as the ratio 
of the fluctuations of some mass tracer, here the AGN (<5agn)> to 
those of the underlying DM mass (5dm), 



b = 



3 AGN 



(10) 



Since the correlation function is generally defined as £(r) = 
(<5(x)<S(x+r)} (the so-called Poisson-process definition), an equiv- 
alent definition of the bias parameter is the square root of the ratio 



Figure 7. Spatial clustering length, rn c > estimated from Eq. (6), as a func- 
tion of X-ray luminosity for the different surveys: The color code and sym- 
bols are the same as in Fig. [4] The black line is the best linear fit to the data, 
excluding the CDF-S (red circles) results. 



of the two-point correlation function of AGN to that of the under- 
lying mass, 



£agn 



1/2 



(ID 



Another related definition, which is the one that will be used in 
the current work, is as the ratio of the variances of the AGN and 
underlying mass density fields, smoothed at some linear scale, 

0"8,AGN 



b = 



C"8,DM 



(12) 



where cts.agn is the rms fluctuations of the AGN density distri- 
bution within spheres of a co-moving radius of 8 /i -1 Mpc, give n 
under the assumption of power-law correlations, bv lPeebles (1980h . 



f"8,AGN = J2(7) 2 (-§- 



(13) 



and as, dm is the variance of the DM density fluctuation field which 
evolves according to (e.g., Peebles 1993): 



08,Dm(z) = 0"8 



D(0) 



(14) 



with D(z) the linear growth factor, which for the concordance A 
cosmology is given by (Peebles 1993) 



D(z) = 



5Q m E{z) 



i: 



(1 + 2/) 



dy 



(15) 



2 J, £ 3 (y) 

Combining Eqs. J 1 2b . dl3t . d 14b . we obtain the cosmological evolu- 
tion of biasing as a function of the power-law clustering parameters, 

'(T 8 D(z 



*>-(?)' 



/ - 



D(0) 



(16) 



Finally, an alternative approach to estimate the bias, free of the 
power-law clustering assumption, is provided by the square-root of 
the ratio of the projected AGN correlation function to the two-halo 
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Table 4. The bias corresponding to the clustering of the Joint sample (1466 
sources) and of the luminosity divided subsamples for each individual X- 
ray survey studied. The corresponding values of the halo mass, in units 
of log 10 (M\ 1 /[h~ 1 Mq]), are estimated using a WMAP7 cosmology, i.e., 
as = 0.81 and Q m = 0.273 and are based on the BPR (4th column) and 
TRK (5th column) bias evolution model. 



l°gl 



b(z) 



BPR 



TRK 



< 42.2 
> 42.2 



41.84 
43.18 



JOINT 
0.98 2.26 ±0.16 

CDF-N 
0.84 1.52 ±0.17 
0.96 2.15 ±0.24 



13.061°;^ 



19 qo+0.29 

,+0.20 
-0.24 



12.S 



jo in + 0.07 



1 2 26 +0,23 

12 7Q+ - 19 
lz -' a -0.21 







CDF-S 






42.5 


42.19 


0.95 4.22 ± 0.47 


14.071^ 


13.86±°- 


> 42.5 


43.45 


1.02 3.75 ± 0.56 


13.80to.20 


13 52 +ai7 






AEGIS 






43 
> 43 


42.75 
43.86 


0.89 1.74 ±0.41 
0.89 2.14 ±0.55 


I2.65i°;« 

13 09+ ' 41 
1J - ua -0 .68 


12 4Q+ - 39 
12.87lE5; 3 5 5 7 






COSMOS 






< 43.5 
> 43.5 


43.30 
44.21 


1.07 2.72 ± 0.50 
1.07 4.12 ±0.96 


i3.24lg;S 
13.86l°:» 


13 04+ ' 24 
i,5 ' u *-0 .33 
13 5Q+0 24 

id.oy_ 34 






eCDF-S 






< 43.2 
> 43.2 


42.86 
43.97 


1.03 2.25 ± 0.83 
1.03 3.00 ±0.97 


12.96+?^ 
13 46+ ' 42 


12 79+ - 47 

lz - ,a -0.97 
1 ^ OO+0.37 
io.zz_ g2 



term of the theoretical projected correlation function, based on the 
Fourier transform of the linear CDM power-spectrum (eg., Allevato 
et al. 2011; Krumpe et al. 2012). These authors have used both of 
the latter two apporaches to estimate the bias and have shown that 
they provide mutually consistent results. 

The CDM structure formation scenario predicts that the 
bias factor is determined by the mass of t he DM halo withi n 
which the extragalactic mass tracer forms (e.g. jMo & White 1996t ). 
In order to assign a characteristic DM halo mass to the es- 
timated bias facto rs we use two bias evolution models, the 
iTinker et al. (2010h (hereafter TRK) and the lBasilakos et al. (20081) 
model ( hereafter BPR). The f ormer, an improvement of the 
original ISheth & Tormen (19991) model, belongs to the so- 
called galaxy merging bias family, which is based on the 
Press & Schechter (19741) formalism, the peak-background split 
i Bardeen et al. 198fj|) and the spherical collapse model (e.g., 
ISheth & Tormen 199dValageas 200°jValageas 201 1]) . The latter is 
an extension of the so-called galaxy conserving bias models, and 
uses linear perturbation theory, and the Friedmann-Lemaitre solu- 
tions to derive a second-order differential equation for the evolution 
of bias, assuming that the tracer and the underlying mass share the 
same dynamics. Furthermore, the model includes the contribution 
of an evolving DM halo population, due to processes like merging. 
Details of these models, as well as a co mparison among five curre nt 
bias evolution models can be found in Pa pageorgiou et al. ("2012b . 



5.1 X-ray AGN Bias 

When applying the above analysis to the clustering results of our 
Joint X-ray sample (with a median redshift (z) ~ 0.98), we find a 
bias factor of b ~ 2.26±0.16. This value corresponds to a DM halo 



Table 5. The values of the X-ray AGN bias factor in redshift bins. Above 
2 = 2 the correlation function was extremely noisy and impossible to de- 
termine. The halo mass shown are estimated using a WMAP7 cosmology, 
i.e., as = 0.81 and f2 m = 0.273 and for three different models of bias 
evolution (see also Papageorgiou et al. 2012). 



# 



2-range 



b(z) logioCAfh/fc^Mo) 
BPR TRK 



353 0.00 - 0.67 

354 0.67 - 0.96 
352 0.96 - 1.38 
200 1.38 -2.00 



0.488 
0.787 
1.123 
1.613 



1.71 ±0.23 
2.46 ±0.25 
2.46 ±0.23 
3.85 ± 1.55 



1007+O.27 ionn+0.21 
10.01 () 36 io.uu_ 2g 

H0.16 1 o 99+O.14 
IO.ZZ n , a 



13.52_ ig u."_ 0il6 
13.00±g"i| 12.90t°-i° 
13.09+?; 4 ! 



0.22 
13.00l°;« 



mass of log 10 (M h //i- 1 M Q ) ~ 13.06(±0.13) and 13.19(±0.08) 
for the BPR and TRK models, respectively (the first raw of Table 
[4^. We also extend the same analysis to the low-L^ and high-L^ 
AGN subsamples in each X-ray field separately. As expected by 
the dependence of the correlation function on luminosity that we 
presented in Section 4.3, the high-Z^ subsamples provide larger 
bias factors and correspondingly larger DM halo masses (for both 
the BPR and TRK bias models) than those of the low-L^ sub- 
samples. An exception is the CDF-S field, which as we have al- 
rea dy discussed, is h ampered by the presence of a large superclus- 
ten lGilli et al. (2003T ). Thus, the data seem to suggest that high-L^ 
sources inhabit more massive DM halos than low-L x sources. The 
results listed in Table|4]indicate that for an average luminosity dif- 
ference of (<51og 10 L x ) ~ 1 between the high and low-L^ AGN 
subsamples, the corresponding average DM halo mass of the for- 
mer sample is a factor of ~ 3 larger than t hat of the latter. 

Similar results have been found by ICappelluti et al. (2010l) . 
for a sample of X-ray selected AGN from the SWIFT/BAT sur- 
vey. The authors derive for AGN with L x < 10 44 erg s _1 a DM 

10 12 M o , while for 



host halo mass of log 10 (Mh/[/i 1 Mq\ 

AGN with L x > 10 44 erg s" 1 they fi nd a DM host halo mass of 
~ 10 14 /i^Mm. lKrumpe et al. (20 id) also find evidence for lumi- 
nosity dependence in a sample of the ROSAT- All-sky-survey AGN. 
In particular, they find log lo (Af h /[/i _1 M ]) = 11.83 and 13.10 
for their low and high X-ray luminosity AGN sources (below and 



above L, 



0.1-2.4kcV 



_ iq 44 -3 er g g 1^ respectively. 



5.2 Bias Evolution 

In order to investigate the redshift evolution of the bias factor we 
split the whole sample in four redshift bins. The resulting val- 
ues of the bias factor in each bin, as well the typical DM halo 
masses, based on the two previously mentioned models, are shown 
in Table [5] Fig. [8] (left panel) shows our determination of the bias 
in the different redshift bins together with the c orresponding re- 
sults of the XMM-CO SMOS dAilevato et al. 201 ll). the XBO OTES 
dStarikova et al. 201 lb. the ROSAT dKrumpe et al. 2010|) and 
SWIFT/BAT dCappelluti et al. (2010h ) surveys. The increase of the 
bias factor with increasing redshift is evident, as well as the con- 
sistency of the derived bias values of the different studies at the 
corr esponding redshifts. 

I Allevato et al. (201lb has argued that models of bias evolution 
based on a single host halo mass cannot explain their results, but 
rather there are indications for different halo masses hosting X- 
ray AGN at the different redshifts. We check their suggestion by 
applying the previously mentioned bias-evolution models, and at- 
tempting to fit a single halo mass jointly to all the previously men- 
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tioned bias data. Indeed no such model can fit adequately well the 
results (the resulting reduced x 2 is J> 1.6). The cause of this failure 
is mostly the deviant bias values around z ~ 0.75. If this redshift 
bin is excluded, then the bias data can indeed be fitted by a sin- 
gle host halo mass with log 10 (A/ h /[/i _1 M Q ]) ~ 13.06 ± 0.04, 
X 2 /df= 0.24 and a present epoch bias of 6(0) = 1.1 (black con- 
tinuous line in Fig. [8} for the BPR model. For the TRK model 
we find log 10 (M h /[/i _1 M Q ]) ~ 12.94 ± 0.05, x 2 /df= 0.48 
and 6(0) = 1.18 (black dashed line in Fig. [8). If we use how- 
ever, only the three deviant bias values around z ~ 0.75, we 
find log 10 (Mh/[/i _1 M Q ]) ~ 13.56 ± 0.09, x 2 /df= 0.29 and 
6(0) = 1.24 (magenta continuous line) for the BPR model, and 
Iog 10 (M h /[/i _1 M o ]) =± 13.28 ± 0.10, x 2 /df= 0.29 and 6(0) = 
1.39 (magenta dashed line) for the TRK model. It is interesting to 
point out that the bias values of the current study (red pentagons), 
when excluding the deviant point, provide exactly the same fit as 
the one obtained previously, using all the determinations of the 
bias evolution values (again excluding the deviant points). How- 
ever, the uncertainty is, as expected, larger (for example, for the 
BPR model we obtain log lo (Mh/[ft _1 M0]) ~ 13.08±0.14 with 
a x 2 /df= 0.48). 



5.3 X-ray versus Optically selected AGN 

We also compare in the right panel of Fig.[8]our X-ray AGN results 
with the bias of optically selecte d QSO, based on the 2dF spec- 
troscopic sample of 20000 QSOs dCroom et al. 20051) . on an SDSS 
(DR5) QSO sam ple of ~30000 spectroscopic QSO with z< 2.2 
dRoss et al. 20091) and on a homogeneous sample of ~38000 SDSS 
QSOs with 0.1 < z 5 l lShen et al. 20091) . Since the X-ray data do 
not extend to such high redshifts and in order to consistently com- 
pare optical and X-ray bias data, we choose to exclude the very 
high-z (3< z < 4.5) optical bias data of Shen et al. Note also that 
since each of the optical QSO studies have used a slightly different 
flat ACDM background cosmology, the bias values shown in Fig. [8] 
have been scaled to the WMAP7 cos mology, following the pre- 
scription of Papageor giou et al. (20121) . The results of each of the 
different studies are indicated in the Figure with distinct colours 
and symbol types. 

The BPR model that fits simultaneously all the optical QSO 
bias data, provides a DM halo mass of log lo (A/h/[ft _1 M0]) ~ 
12.50 ± 0.05 with x 2 /df= 0.68. The corresponding values for the 
TRK bias model are practically the same: \og 10 (Mi l /[h~ 1 Mq]) ~ 
12.49 ± 0.04 with x 2 /df= 0.64. It is therefore evident that 
the X-ray selected AGN bias values, and the corresponding DM 
halo masses, are significantly larger than those of the optical 
QSO, indicating that they probably constitute separate families 
of AGN, each proba bly having a distinct BH fueling mechanism. 
lAllevato et al. (20l"ll) also conclude that the DM halos in which X- 
ray selected AGN reside are more massive than those of optically 
selected QSOs. Further below we will provide a tentative expla- 
nation for such a difference between the host halos of optical and 
X-ray selected AGN. 



6 COMPARISON WITH AGN ACCRETION MODELS 

Our analysis provides new insights into the host DM halos of mod- 
erate luminosity X-ray selected AGN. By fitting two different bias 
evolution models, we find that the AGN in our samples inhabit DM 
halos with masses of \og 10 (Mi 1 /[h~ 1 Mq]) ~ 13.1 ± 0.1. The es- 
timated halo mass is significantly higher than the typical halo mass 



of optically selected QSO [log 10 (Mh/[/i _1 M o ]) ~ 12.5; 
see Fig. [8J. Similar res ults are also found in other X - 
ray studies (e.g. ICoil et al. 2009lStarikova et al. 201ll) . 
dAllevato et al. 201 llMountrichas & Georgakakis 20121) and 
therefore the DM halos masses estimated in most X-ray surveys 
are consistent mostly with those of elliptical/red galaxies. This 
implies that the fueling mechanism of the moderate luminosity 
X-ray selected AGN and the more luminous optically selected 
QSOs is different. 

Semi-analytic models that assume major mergers as the 
main mechanism for triggering AGN activity predict that the 
mass of DM halos that host AGN is lower than that estimated 
from observations of moderate lu minosity X-ray AGN (e.g., 
iMarulli et al. 2009lBonoli et al. 20091) . The theoretical predictions 
from these models are more consistent with the halo mass of 
optically selected QSOs. Additional evidence against the major 
mergers scenario in moderate luminosity X-ray AGN comes 
from the morphol ogical analysis of AGN in the AEGIS an d 
COSMOS surveys dGeorgakakis et al. 2009lSilverman et al. 20091) . 
These authors find that X-ray selected AGN span a large range of 
environments and morphologies with roughly equal numbers of 
bulges, spirals and morphologically disturbed galaxies (see a lso 
ISchawinski et al. 201llKocevski et al. 2012lCisternas et al. 201 lh. 
Also stochastic accretion models dHopkins & Hernquist 20061) 
cannot explain large DM halos for moderate luminosity AGN. 
According to these models, disc instabilities and minor interactions 
feed at high accretion rates relatively small black holes. These 
models predict that moderate luminosity AGN reside in low density 
environments, similar to those of blue star-forming galaxies. 

Regarding the luminosity dependence of the cluster- 
ing of X-ray selected AGN (and consequent luminosity de- 
pendence of the host halo mass) various theoretical mod- 
els of black hole and galaxy co-evolution predict only a 
weak dependence of AGN clustering on luminosity, (e .g., 
iLidz et al. 200dHopkins et al. 200 fepkins & Hernquist 20061) . In 
these models, both the bright and faint AGN reside in similar mass 
halos. The bright AGN correspond to black holes that radiate at 
their peak luminosities, i.e. at accretion rates close to the Edding- 
ton ones. In contrast, the faint end of the AGN luminosity function 
consists of AGN in dimmer (late) phases in their evolution. 

Our d erived DM halo masses of X-ray se lected AGN could 
suggest, as iMountrichas & Georgakakis (20121) also point out, a 
fueling mechanism similar to that of Ciotti & Ostriker (200 lh . 
The authors suggest that stellar winds in early-type galaxies 
provide the gas supply to the black hole (see also Kauffmann 
& Hec kman 2009). Alte rnatively, another plausible model is 
that of iBower et al. (20061) , where gas is accreted onto BHs di- 
rectly from the hot halo of the galaxy. Fanid akis et al. (20121) . re- 
cently presented a calculation for accreting black holes within the 
IBower et al. (20061) semi-analytic model, where AGN activity is 
coupled to the evolution of the host galaxy. The authors assume 
that black holes grow via accretion triggered by galaxy mergers 
and disk instabilities (starburst mode), as well as accretion of hot 
gas from the halo of the galaxy (hot-halo mode). In Fig(9] we com- 
pare our halo mass estimates for the individual subsam ples in Ta- 
ble 5 with the predictions of the Fanidakis et al. (20121) model for 
the L x — Mh correlation at z = 0.5 — 1.5. Their calculations as- 
sume a flat cosmology with Q m — 0.25, ilb = 0.045, erg = 0.9 
and h = 0.7, and for consistency, we re-estimated the DM halo 
mass values according to their adopted cosmology. 

As illustrated by the left panel of Fig. [9] the [Fanidakis et all 
model predicts a strong X-ray luminosity - halo mass correla- 
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Figure 8. Left Panel: Comparison of the redshift evolution of the bias parameter of our Joint X-ray sample (red pentagons), with that of other X-ray surveys. 
The point and color coding of the results based on different surveys is indicated by the corresponding labels. The theoretical bias evolution models, fitted to 
the data for a single DM halo mass, are shown as continuous or dashed lines: The black continuous line is the best fit BPR bias evolution model to all the 
bias data, excluding the deviant points at z ~ 0.75, while the dashed line is the corresponding best fit TRK model. The magenta continuous and dashed lines 
are the fits of the corresponding two models to the bias values around z ~ 0.75. In the inset panel we show the best fit DM halo mass value as a function 
of x 2 — Xmin f° r tne BPR model. A value of x 2 — Xmin = 1 or 4 correspond to the 1 and 2a uncertainty in the fitted parameter (for the case of one free 
parameter). Right Panel: Comparison of the redshift evolution of the bias parameter of the Joint X-ray sample (red pentagons), to that of optical QSOs. The 
point and color coding of the results based on different surveys is indicated by the corresponding labels. The curves correspond to the BPR bias evolution 
model fitted to the X-ray data (red continuous line) and to the optical QSO (black dashed line). Note that the plotted bias values are rescaled to the WMAP7 
cosmology, as in Papageorgiou et al. (2012). 
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Figure 9. The median X-ray luminosity (0.5-8 keV) of the luminosity separated subsamples (from Table 4) v s. the DM halo mass p redicted by the TRK 
bias model. The color coding is same as in Fig. [4] For comparison we show the theoretical predictions of the iFanidakis et al. (20121) black hole model at 
z = 0.5 — 1.5. Solid lines indicate the median of the log L x — log M h correlation for AGN accreting from the hot halo (hot-halo mode; left panel), or during 
disc instabilities and mergers (starburst mode; middle panel) and all AGN in their model (right panel). The dashed lines indicate the 10 — 90 percentiles around 
the median. 
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tion for AGN accreting in the hot-halo mode. In contrast, the 
halo mass of AGN in the starburst mode shows a very modest in- 
crease to X-ray luminosities of Lx — 10 44 ergs~ 4 , beyond which 
it remains constant and approximately equal to 10 12 Mq (middle 
panel). When considering all the AGN in the model, the shape 
of the correlation is dominated by the hot-halo mode in the low- 
luminosity and by the starburst mode in the high-luminosity regime 
(right panel). Thus, the resulting correlation increases steeply in 
the moderate luminosity regime and then sharply declines and flat- 
tens for the brightest luminosities. Our X-ray selected AGN results 
appear to be in very good agreement with the predictions of the 
Fanidakis et al. model. The data lie in the region of the Lx — Mh 
plane, which is dominated by the hot-halo mode. This suggests that 
intermediate luminosity AGN are preferentially powered by accre- 
tion of gas from the hot halo. Their model further proposes that 
sources with Lx > 10 44 ergs _1 , expected to be visible as QSOs, 
live in haloes with masses of- 1O 12 M as suggested by clustering 
studies of optically selected QSOs (see Fig. 8, right panel). 

The distinctive environmental de pendence of the starburst and 
hot-halo mode in the lFanidakis et al J model is linked to the cooling 
properties of the gas in DM haloes l lWhite & Frenk 199ll) . In the 
starburst mode, disc instabilities and galaxy mergers take place in 
gas rich environments. In the standard ACDM paradigm, this cor- 
responds to intermediate mass DM haloes, where gas cools rapidly. 
In more massive haloes, the gas is in quasi-hydrostatic equilib- 
rium and therefore characterized by a significant lower cooling ef- 
ficiency. In addition, in this mass regime AGN feedback reheats 
the gas in the halo, shutting off the cooling completely. AGN ac- 
tivity maximizes at the transition between the rapid-cooling and 
quasi-hydrostatic regime, which in the Fanidakis et al. model is 
found to take place at ~ 1O 12 M0. In haloes more massive 
than ~ 1O 12 M0, black holes accrete low-density gas directly from 
the hot halo. In this case, the authors assume that the amount of 
gas that is accreted onto the BH is determined by, and in fact 
is p roportional to the cooling luminosity of the gas in the halo 
(see lBower et al. 200dFanidakis et al. 201 ll) . Since the cooling lu- 
minosity increases with increasing halo mass, their model predicts 
a strong correlation between accretion luminosity and DM halo 
mass. The good agreement of the predictions of these accretion pre- 
scriptions with our analysis, but also with the clustering results of 
optically selected QSOs, strengthens the idea of two distinct fu- 
eling modes; the hot-halo mode in moderate X-ray AGN and the 
starburst mode in the more luminous optical QSOs. 



7 SUMMARY 

We use 1466 X-ray selected AGN in the 0.5-8 keV band with spec- 
troscopic redshifts spanning the redshift interval < z < 3, with 
a median of z = 0.976. We derive a spatial clustering length of 
r = 7.2 ± 0.6ft" 1 Mpc and a slope of 7 = 1.47 ± 0.12. The 
corresponding clustering length for the nominal slope of 7 = 1.8 
is ro = 6.5 ± 0.4 h^ 1 Mpc. The above clustering length corre- 
sponds to a bias of b(z) — 2.26 ± 0.16, translating to a mass of the 
host DM halo of M h ~ 1.3 x 10 13 h~ 1 M Q . The derived bias and 
the corresponding host halo mass of X-ray selected AGN is signif- 
icantly higher than that of optically selected AGN. This may imply 
a different fueling mechanism for X-ray selected AGN in compar- 
ison to optically selected AGN. We also find indications for a de- 
pendence of the clustering strength on the X-ray luminosity, when 
we consider AGN at a redshift of z ~ 1. The more X-ray lumi- 
nous the sources are the larger, typically, the clustering length and 



hence the higher the DM host halo mass is. For an average luminos- 
ity difference of (<51og 10 L x ) ~ 1, between the high and low-L^ 
AGN subsamples, the corresponding average DM halo mass of the 
former sample is a factor of ~ 3 larger than that of the latter. 

The luminosity-dependent clustering that we find seems to be 
in favor of a black hole accretion model where moderate luminosity 
X-ray AGN, at 2 ~ 1, are fueled mostly by accretion from the 
hot halo around the host galaxy. The good agreement of the model 
predictions with similar studies in the optical further suggests that 
luminous quasar activity is triggered in lower-mass halo (and gas- 
rich) environments by disk instabilities and galaxy mergers. 
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